Probing the largest scale structure in the universe with polarization map of galaxy 

clusters 
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We introduce a new formalism to describe the polarization signal of galaxy clusters on the whole 
sky. We show that a sparsely sampled, half-sky map of the cluster polarization signal at z ~ 1 would 
allow to better characterize the very large scale density fluctuations. While the horizon length is 
smaller in the past, two other competing effects significantly remove the contribution of the small 
scale fluctuations from the quadrupole polarization pattern at z ~ 1. For the standard ACDM 
universe with vanishing tensor mode, the quadrupole moment of the temperature anisotropy probed 
by WMAP is expected to have a ~ 32% contribution from fluctuations on scales below 6.3/i~ Gpc. 
This percentage would be reduced to ~ 2% level for the quadrupole moment of polarization pattern 
at z ~ 1. A cluster polarization map at z ~ 1 would shed light on the potentially anomalous features 
of the largest scale structure in the observable universe. 



1 ) Introduction One of the most intriguing features of 
the temperature anisotropy measured by WMAP [jj] is 
the extremely low value of the estimated quadrupole mo- 
ment. According to the WMAP team, the probability 
of observing such a low value, given their best-fit cos- 
mological model, is ~ 0.7% .2]. Independent analysis of 
the WMAP data have somewhat reduced the significance 
of this finding (see e.g. Q), but have also evidenced 
other curious results, such as the indication of a pre- 
ferred direction for the quadrupole and octopole modes 
PJ ■ These discoveries have raised speculations on possible 
non-standard topology of the Universe 0, or inflation- 
ary physics (see e.g. and at the same time, evidenced 
the importance of independent observational methods to 
probe the largest scale structure in the Universe. 

The low-Z moments of the present temperature 
anisotropies are usually assumed to probe the largest 
scales, but in fact they carry information about perturba- 
tions on a fairly broad range of scales. It would be highly 
preferable to find an observable which probes the largest 
scales with less contamination from intermediate ones. 
In this Letter we point out that the polarization signal 
from high-redshift clusters could be such observable. 

Light scattered off a free electron is polarized, if the 
electron sees a quadrupole anisotropy of its incident light. 
This effect gives rise to interesting observable phenom- 
ena, like, for example, the large-scale polarization signal 
in a reionized Universe 0, 0. The same effect is re- 
sponsible for the major polarization signal in Sunyaeve- 
Zeldovich (SZ) galaxy clusters. As first pointed out by 
, the polarization pattern of clusters at different red- 
shifts may shed light on the quadrupole amplitude at 
earlier times and at different positions (see also 0] in 
relation to the low observed quadrupole moment). The 
cluster polarization has also been investigated as a po- 
tential tool to determine the nature of dark energy that 
is sensitive to the Integrated Sachs- Wolfe effect 0, Q| . 

In this Letter we develop a transparent formalism for 
the calculation of the relevant observable quantities of 
the cluster's polarization signal, and show that the po- 



larization pattern at z ~ 1 could be a promising probe 
of the largest scale structure in the universe, compared 
with the quadrupole temperature anisotropies observed 
today. This finding contrasts with the naive expectation 
that the polarization signal of distant clusters are gener- 
ated by fluctuations whose spatial scales are well within 
the horizon at the present [12j. We also estimate how 
well the dark energy equation of state can be measured 
by the sole cluster polarization signal. 

2) Formulation In this section we review the formal- 
ism to analyze the polarization signals of galaxy clusters. 
Throughout this Letter, we only discuss linear scalar per- 
turbations in a flat background universe. We first intro- 
duce a fixed spherical coordinate system Co centered on 
an observer at z = 0. In the coordinate system Co wc 
denote the angular variables with Q — (9, (f>) and use red- 
shift z as a radial coordinate for our past light cone. Let 
us consider the linear polarization generated by the lo- 
cal temperature anisotropy seen by a cluster labeled by 
i at a redshift Zi and direction Qi . Such polarization can 
be expressed as a linear combination of the quadrupole 
anisotropy CL2m{^ii z i) seen by the cluster as 
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Xi = ~r~y~r — = ^ a2 m (^i,Zi) 2 Y 2m (^i), (1) 

V ' 771— — 2 

where F(z) = -^/IQT CM b{z) is a normalization fac- 
tor (Tcmb(z): CMB temperature at redshift z), r C i is 
the effective optical depth of the cluster, and iYi m iS^) is 
the spin- weighted spherical harmonics [lj, [l5j defined in 
the Co system. The coefficient 0,2m (^i,Zj) in eq.QJ is 
defined in a spherical coordinate system Cj whose ori- 
entation is obtained by a parallel transport of Co to the 
cluster's position (fij,^). We shall assume that we can 
make a three dimensional polarization map X(Q,z) on 
our past light cone by observing many clusters at differ- 
ent redshifts and directions. In eq.Q we only include 
the primary temperature quadrupole anisotropy as the 
source for the cluster polarization, since the contribu- 
tion of the secondary polarized incident light is expected 
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to be much weaker. We also neglect the effects of the 
peculiar velocity field 0, |n| , assuming that the typical 
comoving distance between the surveyed clusters is larger 
than the correlation length of the peculiar velocity field 
< 50/i -1 Mpc 0- To make the map X(Sl,z) we need 
to estimate the optical depth r ci of each cluster. This 
would be performed by using other observations like the 
SZ spectral distortion or the X-rays. 

One way to characterize the statistical properties of 
the polarization map X(£l, z) is through its correlation 
function (X(Qi,Zi)X(ilj,Zj)). This function is written 
in terms of the correlation (a-> m {Qi, Zj)a?. m i f O, , Zj)) for 
the information given at two positions i and j [l2j. But 
this representation has two disadvantages. First, it is 
given by different frames of reference. Second, it is not a 
diagonal matrix with respect to m and m'. As a result its 
expression is very complicated and hard to relate to basic 
theoretical inputs (e.g. the primordial power spectrum). 

Here, we extensively use the properties of the spin- 
weighted spherical harmonics for analyzing tensor quan- 
tities. This kind of approach is widely used in both CMB 
and weak lensing studies 0,0,0], and especially useful 
for dealing with statistically isotropic fluctuations that 
are analyzed in this Letter. Our goal here is to write the 
map X(Q, z) in the orthonormal form; 

oc I 

X(Q,z)=Y^ J2 2 Y lm(n)bi m (z), (2) 

1=2 m=-l 

and to present the basic formulas that relate the coef- 
ficients bi m (z) to the spectrum P(k) of the primordial 
density fluctuations. In order to cast eq.Q) as in J5J, 
we first separate the radial information z and the direc- 
tional information Q for the coefficient a2 m (0,, z). Then 
we combine the angular information with that of the spin- 
weighted harmonics 2^2m(^) in eq.JQ) to get the expan- 
sion with the orthonormal angular basis 2^m(^) as in 
eq.©. The latter process is similar to the unification 
of the spin and orbital angular momentum in Quantum 
Mechanics (e.g. ^H), whose application to CMB is dis- 
cussed in Hu & White 0] (see also H3)- Following the 
analysis for the CMB polarization in a reionized universe 
[Ullll, we find 

(b lm (z)tf, m ,(z')) = 8 u >5 mm ,{Anf (3) 
x J ^(P(k)k- 1 )h l (k,z)hi{k,z'), 

where P(k) = Ak n (A: a normalization factor) is the 
primordial power spectrum with n = 1 corresponding to 
the scale invariant spectrum. Due to the assumption of 
the statistical isotropy, the correlation (3) has diagonal 
form. This expression is given for E-mode (electronic par- 
ity) polarization generated by scalar perturbations that 
would not produce B-mode (magnetic parity) polariza- 
tion. Tensor (gravitational wave) or vector perturbations 



can generate both E, and B-mode polarizations, and we 
can determine or constrain their amplitudes by measure- 
ing B-mode polarization. 

In eq.© the function hi(k, z) is defined as 

h l (k,z)=A 2 (z,k)f l [k( T (0)-T(z))], (4) 

where r is the conformal time, A2(z,k) is the transfer 
function for the quadrupole temperature anisotropies at 
a given z, and /;[/c(t(0) — t(z))} is a geometrical factor 
which relates the scale of the fluctuation with the dis- 
tance of the cluster from the observer. We will return to 
this factor later on. 

The transfer function A 2 (z,k) is related to the linear 
growth rate D and to the scale factor a as: 

3 

A 2 (z,k) = —j 2 [k(T(z) ~ T(z rec ))} 

+ §£^ W ^)-r)]£(5«)(» 

where j2{x) is the spherical Bessel function and r rec is 
the conformal time at recombination z rec ~ 1100. Note 
that the function A 2 (z,fc) represents the evolution and 
the projection effects of each Fourier mode k, but it does 
not depend on the power spectrum. The first term on 
the r.h.s. of eq. JSJ) is the Sachs- Wolfe (SW) effect, and 
conveys information including the largest scale structure 
(comparable to the horizon size) at redshift z. The sec- 
ond term is the Integrated Sachs- Wolfe (ISW) effect, and 
is sensitive to the recent expansion history of the uni- 
verse. The ISW term typically probes smaller scales than 
the horizon size. 

The function fi(x) in eq.Q represents projection ef- 
fects for scales of the order of the cluster's distance from 
the observer and it is expressed in terms of spherical 

Bessel functions as f t (x) = \J ^-2)\ 30^- 

Note that eq.JQ) only contains the expression of the 
local quadrupole (/ = 2) mode at the cluster's location, 
while the expression in eq. J3J) also contains higher modes 
(I > 3). This is due to the power transfer from I = 2 to 
I > 3 which is caused by the spin-orbit angular momen- 
tum coupling. The function fi(x) regulates such transfer, 
preserving the total power (J^ 2 fi( x ) 2 (2l + 1) = 1). 

We can now build a natural estimator of the total an- 
gular power spectrum for each i-mode from the observed 
map as follows: Hi(z) = Y^ m =~i \hm{z)\ 2 . We assume 
that the primordial potential fluctuations are random 
Gaussian distributed. Then, it is straightforward to cal- 
culate the covariance of the spectrum 

Coo{Hi{z),H v {z!)) = 2<W(2l+l) (T (M*)CCO>) 

(6) 

where ((bi m (z)b* m (z'))) is given by eq.©. At z = z' 
this expression trivially allows to evaluate the cosmic 
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variance for Hi(z), leading to the familiar expression 
Cov(H l (z) 1 H l (z) = 2/{2l + l)H l (z) 2 . 

3) Results In this section we present the general fea- 
tures of the power spectrum Hi(z). We first discuss how 
the observed quadrupole spectrum H 2 (z) is related to 
the matter fluctuations at different wave number k. As 
shown in eq.® the term hiik^z) 1 represents the weight 
for the function P(k) / k which is constant for the scale in- 
variant model. Therefore, we shall use h2(k, z) as a mea- 
sure of the spatial scale probed by the spectrum H 2 (z). 

In figure 1 we plot the function h 2 (k,z) at three dif- 
ferent epochs z = 0, 0.5 and 1 for our fiducial cosmolog- 
ical parameters Q m = 0.3 and fl\ = 0.7. Note that the 
curve for z — can also be a regarded as the weight for 
the quadrupole moment of the temperature anisotropics 
observed today. This curve shows that the temperature 
quadrupole receives a contribution from fluctuations on a 
broad range of scales ~ 10 _4 /iMpc _1 to ~ I0~ 2 /iMpc -1 . 
The wiggles in the curve reflect the oscillatory behavior of 
the spherical Bessel function in the first term on the r.h.s. 
of eq.JSJ) (related to the SW effect), while the negative 
mean value around 10 _3 /iMpc _1 < k < 10~ 2 /iMpc _1 is 
due to the second term (ISW effect). 

Let us discuss the redshift evolution of the function 
h 2 (k, z). The wave number of the first peak of the curves 
h 2 (k,z) is determined by the horizon scale at each red- 
shift. As expected, such wave number increases as we 
move to higher redshift. This scale, however, does not 
change much at z < 1: for our fiducial cosmological 
model the comoving length of the horizon is reduced only 
by 14% (24%) at z = 0.5 (z = 1) compared with the hori- 
zon size at z = 0. 

Let us now discuss two other important effects that 
tend to increase the weight of large spatial scales probed 
by H2{z) at redshift z ~ 1. The first effect is due to 
the presence of dark energy. In general, our Universe 
becomes very close to the Einstein de-Sitter one at a rel- 
atively low redshift z ~ 1. As a consequence, the weight 
associated with scales smaller than the horizon (which 
are the ones typically affected by the ISW effect) is sig- 
nificantly reduced at z ~ 1. 

The other important effect that suppresses h 2 (k,z = 1) 
at small scale is the transfer of power from I = 2 mode to 
higher modes. This effect is characterized by the func- 
tion f2(x) that has a profile f 2 (x) ~ 1 at i < 1 and 
f 2 (x) ~ — 15sin(cc)/a; 3 at x > 2. The typical scale be- 
low which this suppression occurs is proportional to the 
distance to the cluster. However, as we are dealing with 
our past light cone, this distance coincides with the dif- 
ference in the horizon sizes at the present and at the 
cluster's redshift. As noted above, this length is 24 % of 
the present horizon for clusters at z — 1, so that fluctua- 
tions on scales smaller than this one cannot make a im- 
portant contribution to the observed quadrupole moment 
H2(z = 1). The present result is straightforward in our 
new formalism, while it has been overlooked in previous 
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FIG. 1: Redshift evolution of the function hiik,z) for the 
fiducial cosmological model with Q m = 0.3 and Qx = 0.7. 
Overall scale is irrelevant here. There are three effects that 
change the contribution of fc-mode to the quadrupole pattern 
H2(z). Arrows show how these effects change the shape of 
the curve h,2(k, z) with increasing z. 



works 12j because of the complications introduced by the 
use of multiple coordinate systems (see also 0). From 
figure 1, it is apparent that the polarization map around 
z = 1 would be a powerful tool to study the largest scale 
fluctuations avoiding the small scale contamination. 

We will now attempt a quantitative analyses aimed to 
determine which wave number k provides the dominant 
contribution to the function hi(k,z). We define a func- 
tion R(z) that is the mean of log 10 [fc] weighted by the 

weight h 2 [k,zy as follows R(z) = d r . 

If we define log 10 [fc*] = R(z), the wavelength fc* can 
be regarded as a typical scale probed by the moment 
H 2 (z). We obtained R(0) = -3.16, i?(0.6) = -3.42, 
R(1.0) = -3.41 and R(2.0) = -3.38 (with k is in units 
of ZiMpc -1 ). Roughly speaking, we can reduce the ef- 
fective wave number by a factor of ~ 2 by using a map 
at z ~ 1, compared with the quadrupole temperature 
anisotropics observed today. For the sake of comparison 
we also calculated the same kind of quantity for the tem- 
perature anisotropics at z = 0, and found R = —3.08 
(R = -3.00) for 1 = 3(1 = A). 

In addition to R(z), we also studied the contributions 
of fluctuations below and above k = 10 _3 /iMpc~ 1 to 
the spectrum H2(z). In figure 1 the wave number k — 
10 _3 /iMpc _1 is given by the vertical long-dashed line, 
and corresponds to ~ 6.3/i _1 Gpc in real scale. At z = 
about 32% of the power is coming from k > 10 _3 /iMpc _1 , 
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FIG. 2: Redshift evolution of the spectrum Hi(z) for param- 
eters fi m = 0.3, Q\ = 0.7 and n = 1. 



but its contribution decreases to ~ 2% at z = 1. We 
conclude that the quadrupole of the cluster polarization 
signal at z ~ 1 would allow to probe the large scale power 
of the Universe in a cleaner way. 

So far, we mainly discussed the quadrupole mode. As 
we commented earlier, we expect higher modes (I > 3) 
to be generated by the spin-orbit coupling. In figure 2 
we show the spectrum Hi(z) as a function of redshift. At 
z = the polarization map X simply reflects our local 
quadrupole mode so that -H)>3 = 0. The total power 
Hi( z ) (dashed curve) is the averaged local temper- 
ature quadrupole moment J2 m ( a 2m a 2m) & f each redshift 
(see eq.Q). At z > 1 it becomes a constant value due to 
the scale invariance of the matter power spectrum. The 
quadrupole mode H2{z) decreases continuously with in- 
creasing redshift due to the power transfer to higher-^ 
modes. However, the I = 2 mode remains the largest sig- 
nal at redshifts z < 2 where the cluster polarization map 
would be observationally available. 

We shall now consider if the redshit dependence of the 
polarization signal can be used to constrain cosmology. 
In particular we investigate possible constraints on the 
equation of state of dark energy. We evaluate the pa- 
rameter estimation errors using the Fisher matrix ap- 
proach applied to H%{z). We combine information on 
H2(z) from different redshift shells between < z < 2 
binned in Sz = 0.2 intervals, and assume cosmic variance 
with appropriate bin correlations (see eq.®, and also 
[l^)- as the sole source of the error. 

We constrain a single fitting parameter w = 
P/p ^constant around the fiducial model (f2 m = 0.3, 
£l\ = 0.7 and n — 1), and find Aw ~ 0.6 (lc). Such 



a large error is due to the large cosmic variance in each 
redshift bin and to the strong correlations between bins. 
Indeed the SW effect is not important for dark energy 
studies, but it contributes to the above correlation. Thus, 
in the attempt of removing the SW contribution, we 
also applied the Fisher matrix approach to the estima- 
tor E m (\b2m(z) - ^ m (0)| 2 ), and found Aw ~ 0.4 (la). 
In order to improve this result we could also include, in 
principle, information of the higher order modes / > 3. In 
an actual observational analysis, we need a lot of efforts 
to detect these weak signals around at z < 1 where the 
effect of the dark energy is important (see figure 2). We 
conclude that cluster polarization alone would not be a 
powerful observational method to constrain dark energy. 

Finally let us make an order of magnitude estimate 
for observational requirement needed to measure the mo- 
ment H 2 (z) at z ~ 1 (see also 0). Suppose we measure 
the polarization signal for N clusters distributed on the 
whole sky at z ~ 1 with an observational error a on the 
polarization intensity for each cluster. The total error 
for Hi is AHi = y/2(2l + l)[Ht/(2l + 1) + 4ttct 2 /A] with 
the first term representing the cosmic variance and the 
second term being the measurement error. For I = 2 the 
two errors are equal if a ~ 0.3(t c /10- 2 )(A/10 3 ) 1 /VK. A 
part from the sensitivity considerations, there is also the 
issue of separating the cluster polarization signal from 
the other competing ones, like, for example, the CMB 
lensing. This is best achieved if the cluster is spatially 
resolved, which typically requires observations with res- 
olution of ~ 1 arcmin. Planck, for example, is an all sky 
survey mission with sensitivity ~ 5fiK per pixel (for com- 
bined polarization channels and 14 months integration) 
and ~ 5 arcmin angular resolution. Because of the low 
resolution, it may not be the most suited experiment to 
detect H2 around z ~ 1. 

As the polarization pattern is dominated by the low-^ 
modes, a fine sky sampling is not necessary to map it. 
Furthermore, in order to probe only / = 2 mode at z < 1 
where I > 3 modes arc weak, wc just need to observe half 
of the sky (27r[sr]) due to the parity symmetry of scalar 
perturbations (though, in this case, the power transfer 
in figure 1 does not work). A ground-based telescope 
performing targeted cluster observations in small areas 
sparsely distributed on half sky may be an adequate tool 
to achieve this task. 
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